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Q ■ Abstract 

C^, We analyse the structure of current-current correlators in coordinate space in large 

Q^l Nc limit when the corresponding spectral density takes the form of an infinite sum 

^ ■ over hadron poles. The latter are computed in the QCD string model with quarks at 
the ends, including the lowest states, for all channels. The corresponding correlators 
demonstrate reasonable qualitative agreement with the lattice data without any ad- 

/\^ • ditional fits. Different issues concerning the structure of the short distance OPE are 

j^ . discussed. 



1 Introduction 

Correlators of hadron currents (denoted as CC in what follows) 

P(^)(a;) = (j(^)(x)j(^)(0)t) (1) 

are the basic elements in QCD which have been studied in the framework of perturbation 
theory [1^ and also using the Operator Product Expansion (OPE) [2^. More recently CC 
have been the subject of exciting lattice calculations |H1I1] as a tool to study nonperturbative 
vacuum configurations. Since then it has been realized that CC can provide information of 
two sorts: a) one is encoded in the hadron mass poles m^ and quark coupling constants c„ 
and can be associated with the long distance dynamics (LDD), b) another kind of information 
can be extracted from the behavior of V^^\x) at small x (or at large Q^ from the Fourier 
transform V^'^\Q)) with the help of the OPE (see, e.g. [5]) where the coefficients represent 
vacuum matrix elements of local field operators. Since OPE is assumed to be valid at small 
distances, one can associate it with what we shall call the short distance dynamics (SDD). 

There is another way to stress these two different facets of CC analysis. In the standard 
Shifman-Vainshtein-Zakharov (SVZ) framework j2] one proceeds in two steps. First, the CC 
is expanded in local products of fields and, second, the matrix elements of the latter over 
(nonperturbative) vacuum states are taken. However, in principle one can go the other way: 
first, to take an average of the CC over the vacuum and, as a second step, to expand the 
result in power series at small distances. In the former case the answer is written in terms of 
the SVZ vacuum condensates (like {G^ ) or (ipip)), while in the latter case one has to deal 
with the spectral parameters (like string tension a, masses etc.). It is often assumed in the 
literature that the results of the two procedures outlined above must coincide identically. 
This is an important element of the SVZ sum rules: keeping only a few lowest terms in 
expansions of both kinds, one puts them equal to each other and with the help of some 
additional mathematical tricks (Borel transformation) the spectral parameters of low lying 
resonances can be expressed in terms of universal vacuum characteristics like quark and 
gluon condensates. From the quahtative point of view, the approach forbids the so called 
}? IQ^ term (see discussion in [B]) since no local gauge-invariant dimension 2 operator exists 
in QCD. On the other hand, it is easy to get nontrivial term of this kind (for example, 
A^ ~ cr) from the spectral expansion point of view. 

The two approaches just described have an important physical difference. The standard 
OPE framework encodes the knowledge about confinement properties of the theory in very 
indirect way, via condensates etc. However, nonzero vacuum expectation value of some local 
gauge-invariant operator by itself does not imply confinement since the latter is essentially 
nonlocal phenomenon linked with the structure of asymptotic states of the theory.^ On 
the other hand, the spectral expansion manifestly takes this aspect into account. One can 



^This is also true on the lattice where confinement criteria such as area law or linearly rising potential 
are physically meaningful only for times/distances larger than some typical nonperturbative scale. 



say that the standard condensate series is surely sensitive to nonperturbative physics of the 
theory but at the same time is " confinement-bhnd" at the level of a few lowest terms, since 
confinement is expected to arise after summation of the whole series. 

In our opinion, there are no solid reasons for the mentioned coincidence between SDD 
and LDD to take place beyond perturbation theory. The key point here is divergence of 
the series in question. This precludes to construct one-to-one correspondence between terms 
of the small distance expansions based on LDD and on SDD. In general case each term of 
the expansion based on LDD is a result of some SDD subseries summation and vice versa 
(see examples in |2j). Moreover, strictly speaking one needs an additional prescription to 
assign physical meaning to the (divergent) sum of matrix elements of local operators. We 
do not address general issues of the OPE structure here, leaving this problem for future 
publications (see also jHj and references therein in this respect). Instead we shall compare 
below the explicit expression of CC in LDD at large N,, with the lattice data [Sill]- These 
data provide an essential information on CC in the range < x < 1.5 fm where both LDD 
and SDD are present and therefore one may check in principle i) the region of validity of 
the standard OPE, ii) the transition region from SDD to LDD and the quark-hadron duality 
(QHD), iii) the validity of specific vacuum models and finally iv) the validity of LDD. 

In the latter case one has large N^ limit where there is a solid knowledge of the spectrum 
properties [H]. Some aspects of OPE at large N^. have been discussed recently in [TD]. Let us 
remind that the masses and residues (properly normalized) of CC, rrin and c„, are stable at 
large Nc limit, and recent lattice data (see ^I] and references therein) confirm that typically 
corrections to stable quantities for iVc = 4, .., 8 are small. ^ Moreover at large Nc one can in 
principle compute the hadron spectrum (i.e. the set of m„, c„) quasiclassically in analj^ic 
form for all n and thus calculate CC explicitly. This was done for the masses (see ^21 and 
references therein) and the residues {^^, see also ^3]) and compared in ^H] with OPE, 
SVZ sum rules and experiment. In doing so one encounters and solves several important 
problems, which we briefly comment on below. 

First, the knowledge of m„, c„, at large Nc allows to find the limit of the CC at large 
Q (small x) and compare it with the quark-partonic expression, thus checking the quark- 
hadron duality which was suggested long ago and studied since then in numerous works (see 
reviews ^H] containing extensive lists of references to the original papers). 

Secondly, the large Nc LDD provides the small x OPE with coefficients depending only 
on the string tension a, and this LDD OPE is to be distinguished from the standard, or 
SDD OPE. These two types of OPE are in principle distinct, however in the vector channel 
both agree with the experimental e'^e' data, as was shown in ^Sj. Many efforts have been 
undertaken to match these OPE's (see, e.g. 11111121120]), but here we do not follow this line. 

Thirdly, one has a unique chance to compare theoretical values of CC in LDD with 
experiment and lattice data and discover its region of validity. 



^In what follows the flavour-singlet channels where corrections can be large, are not considered. 



Since theoretical CC in LDD contain no adjustable parameters at all, and are computed 
in terms of fixed values of a and a^ only (and final forms of CC contain only two fixed masses, 
expressed through a, a^) this comparison has a form of a fixed prediction. 

The main purpose of this letter is to construct explicit expressions for CC in large N^. 
LDD in the Euclidean 4d x-space and to compare those with lattice results from j3] and 
|3]. This is done in sections 2 and 3. The data from |3] we are going to compare with were 
obtained on 16^ x 24 lattice with the spacing 0.17 Fm and /? = 5.7, while the author of j3] 
used the lattice of the total size 1.5 Fm and the spacing 0.13 Fm at /? = 5.9 Both simulations 
were performed in quenched approximation, however the current fermions were taken into 
account differently. An interested reader is encouraged to consult the original papers H] 
for all technical details concerning the simulations. 

In section 3 we calculate the LDD OPE as the Taylor a;-expansion and compare it to the 
standard OPE and lattice data, defining in this way regions of validity of the corresponding 
expansions. There are a few subtle points of difference between OPE in the coordinate space 
and in the momentum space, stated already in the original papers j21 121] and discussed also 
in j22]. We are not going to repeat old arguments here and refer an interested reader to the 
corresponding literature. The fourth section is devoted to discussion of the results and to a 
short summary of possible perspectives. 



2 The model 

The model of zero-width equidistant resonances j23| 1^ 03] has been analyzed in different 
respects in connection with OPE and large Nc QCD, where the main interest has been 
concentrated on the momentum-space representation fIE\ IT^ . We are following a different 
path, keeping ourselves in coordinate space. 

As it is discussed in the Introduction, in order to get an information about nonperturba- 
tive contents of the current correlators in QCD one usually proceeds along one of two ways. 
The first one starts from the small distances - due to asymptotic freedom it is safe to assume 
that the leading term at small enough distances is given by the perturbation theory. There 
must be power corrections however, and at this moment the standard method of SVZ comes 
into play. The second way goes from the large distances. Because the theory in question 
is by assumption confining, at large distances current-current correlator is dominated by 
exchange of the lightest excitation with the particular quantum numbers. When distance 
is decreasing, the contribution of higher states is turning on and at intermediate distances 
the whole tower of states plays role. We choose the latter way to go in what follows. The 
main object of our interest is the current-current correlator in coordinate space given by ([Q), 
where the current j^^\x) = ^p{x)Tc^{x) is defined for the given channel c. We consider only 
flavor nonsinglet charged currents (of the type uTd) in this paper. The matrix Tc carries 
Lorentz, flavor and color indices (with the latter structure being always Sap), corresponding 



to quantum numbers of the channel c. So we have for vector (p - channel), axial (ai - chan- 
nel), pseudoscalar (vr - channel) and scalar (ao - channel) the following expressions for the 
(charged) currents 

J^;\x) = u^,d; jt\x)=u^,^,d\ J^'\x)=m^,d- j^'\x) = ud (2) 

We adopt the standard normalization of ^ |2] , which is different for neutral currents by a 
factor 1/a/2 from the one used in j^]. 

It is convenient to factor out the free part of the correlator (^ and to define the ratio 

i?(^)(x) as 

' freeV^J 

where by definition for the vector and axial channels the sum over Lorentz indices is always 
taken: 

v^^'^'^Hx) = g^^ r^;f\x) (4) 

The free part V%[.^{x) is given by the following expression in the chiral limit: 



j{v,a,p,s) 



iVc 



^/"er'^(^) = ;;zi-(2,2, -1, -1) (5) 



The important relation R^'^\x) has to obey is given by 

limi?('=)(x) = 1 (6) 

Its physical interpretation depends on what point of view - SDD or LDD - one takes. In 
the former case the fact that all effects of interaction - perturbative and nonperturbative 
- switches off at small distances is what is commonly known as asymptotic freedom. The 
interpretation is more subtle in the latter case - the fact that the whole tower of hadron states 
gives total contribution coinciding with that of the free quarks (at small distances) is one of 
manifestations of the quark-hadron duality. From the point of view of the spectral expansion, 
duality is a highly nontrivial phenomenon since it puts severe bounds on the properties of 
{rrin}, {cn} at n — i> oo. In principle, one should be able to derive, rather than postulate quark- 
hadron duality for each given channel starting, for example, from (relativistic) Hamiltonian 
and computing its spectrum in this channel. This problem is beyond the scope of our paper 
and will be considered elsewhere. Instead we will always assume quark-hadron duality in 
what follows and fix the coefficients {c„} to yield duality. In other words, the condition (jH)) 
holds by definition in our approach. 

Another problem is the perturbative corrections. As a matter of principle, the LDD 
procedure remains the same, but one is to include in the spectral sum hybrid states. The 
fact that diagrams of perturbative corrections in confining vacuum are connected to the 
hybrid excitations was discovered in [2^] and discussed in [12], see also [2Z1I2H1- At small 
distances n-gluon hybrid excitation would correspond to the perturbative correction 0(a") 



from the SDD point of view. It is an interesting open problem to establish this relation in 
explicit way and we leave it beyond the scope of this paper. The lattice data we are mainly 
focused on provide no clear indication of such corrections as well as of any effects caused 
by perturbative running.^ The same reasoning is applicable to anomalous dimensions for 
nonconserved currents we take no care of in what follows. 

We find it instructive to explain our attitude, which in some respects is different from 
that of the cited papers jTHJ UHl 1201 ■ Our basic assumption is that QCD exhibits confinement 
in the form of minimal area law in large Nc limit. Since we always have in mind large Nc, 
the picture of zero- width states is justified and the confinement property guarantees the 
spectrum to be discrete. On the other hand, thanks to the area law one is able to solve 
the corresponding Hamiltonian problem for two masses (quarks) connected by the minimal 
string [TTj and to get equidistant spectrum (at least quasiclassically) , as it was done in [T2]. 
Therefore for each channel we are dealing with the corresponding set of poles {rrin} and 
residues {c„} which in a sense completely characterizes the theory. We get for the imaginary 
part of polarization operator 

1 oo 

-lmU{s) = J2cn-5{s-ml) (7) 

^ n=0 

This is the basic expression we are to work with below. 

The mass spectrum is obtained quasiclassically to have the following form ^Hl ^1 ISHj 

m^ = ml + rn^n (8) 

where the quantity rn^ is defined universally for all channels to be m? = Ana as it was found 
by the quasiclassical analysis (see also jM] and references therein) . Here a is physical string 
tension. The residues for the ground state and the excited states are treated differently. As 
for the latter, they are fixed by the requirement of quark-hadron duality. The lowest state 
residue is chosen in different physically motivated ways for different channels (see below). 
Our approach is phenomenological in this respect, since we have not computed all c^ starting 
from some consistent theoretical scheme (it will be done elsewhere). It is worth stressing 
however, that we do not perform any kind of fitting procedure. In some sense, we have no 
fitting parameters at all in our formulas since all quantities like mo, cr take their physical 
values. We do not try to fit the lattice data, instead, we use them for comparison with the 
results of quasiclassical spectrum of large N^ QCD. 



2.1 Vector channel 

For conserved vector current one can write the following expression in momentum space 

^g(g) = ^/^'^^g (^) exp(^ga;) = {q,q, - q'g,,)U'^^\q') (9) 

•^For example, due to asymptotic freedom R^'^\x) should approach unity logarithmically slow when x — > 0, 
which is not seen on the lattice, as figures 1.1 - 4.1 clearly demonstrate. 
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and the determination of 11*^^^ (g^) for —q^ = Q^ > is of interest. 

It was shown in P^ IT^ IT7j that for a system made of relativistic quarks connected by 
the straight-hne string with tension a one has 

m\ = 2Txa{2nr + L) + ml ; Cn = —7^ ■ 47rcr (10) 



which leads to 



where 



\2 /xr °° /7 



m^ = m n + ttIq ; an = rn = Ana (12) 

Despite we work in the large Nc framework, we have kept the factor Nc in front of the second 
term in (jllj) in order to make contact with the asymptotic expression pH|) . The residue of 
the first (p-meson) pole should in principle be calculated in the same dynamical framework, 
which is used for the spectrum calculation.^ The condition ()10|) provides the quark-hadron 
duality in this channel: 



n(^)(-g^)'^^--^iog(^) (13) 

Numerically we use values ml = M'^ = 0.6 GeV^ A^ = 0.047 GeV^ and m^ = Ana = 2.1 
GeV^, corresponding to cr = 0.17 GeV^. The value of A^ is consistent with the one used 
in 1201 ^^d also the one computed in jHU] . Notice that in the latter case we are not to 
take into account the large 0{as) correction which is not seen in lattice simulations. It is 
interesting that the value of the lowest state residue A^ is different from the asymptotic value 
Cn by less than 15% in our case. 

Since we want to compare our model with the lattice results, the Wick rotation has 
to be performed. See Appendix B where all relevant formulas are collected. The resulting 
expression in coordinate space, corresponding to (fTTj) is given by: 

Ri^){z)=iz^K,{z) + —p,{z) (14) 

where dimensionless distance z = xe^tlq and mass ratio h = m? /ml have been introduced 
and 

e = ir -4-7-^ (15) 



ml Atx"^ 



The universal function Pb{z) is given by 



_ [du ( z^ \ l + exp(-&n)(6-l) 

^^(^) = / ^ '"" VTu - "j (l-exp(-H)^ ^''^ 



^This is also true for higher resonances, however, we choose another way and fix those residues by quark- 
hadron duahty for those channels where c„ have not yet been calculated dynamically. Notice that in some 
channels like the vector one, dynamically calculated c„ indeed provide quark-hadron duality (see \il\). 



As it has already been mentioned, we have three parameters in the expression ()14|) : mo, 
which fixes the overall scale of distance, A^ and b characterizing the spectrum, and all of 
them are fixed by their physical values. The resulting plot is shown on Figs. 1.1 and 1.2, to 
be compared with the lattice data of jH] and |3], respectively. One sees reasonable qualitative 
agreement with the lattice data. In fact, the results from P] and from Pj were obtained for 
different lattice parameters (see the cited papers for details) and strictly speaking they do 
not agree with each other at large distances. However the qualitative agreement takes the 
place, and both sets of data reasonably correspond to our simple expression (lE^ 



Of separate interest is the short distance limit of our results. It is straightforward to 
expand (|T^ near z = 0, the answer is (see Appendix B) 



R^'^z) 



4 



^z'' + ^z''\ogz^ + 



+ 



+ 



1 + 



3-2^ 



Qb-Q-b'^) / + 



36-2-62 
3^28 



z^ log z'^ 



(17) 



The first and the second brackets in the r.h.s. of (|17p correspond to the first and the second 
terms in the r.h.s. of (|14j). respectively. The dots stay for higher order terms. 

Notice an absence of a term proportional to z'^ in (fT7|) . It is easy to see by dimensional 
reasons that such term would correspond to X^/Q"^ contribution in the polarization operator 
(fTT|). Indeed, high momentum asymptotics of (fTT|) reads: 



n('')(-Q2 



127r2 




+ ( Air^Xl - -m' 



Ulr 



Q'- 



-+ 



+ 



1 4^1 . 

—nip. H — m m, 

2 ° 2 



1 + 



12 



-m 



4-^A>^) ■^,+0{Q 



Q' 



;i8) 



so it generally contains l/Q"^ term (and in this sense violate the standard OPE where such 
term is absent). In coordinate space, however, this term being multiplied by Q^ produces 
delta-function like contribution which we systematically disregard in this paper. In this sense 
it is impossible to detect l/Q^ term by simulation of the correlator in coordinate space (it 
was noticed in |S21)- However, one may take a different attitude and put the requirement 
A^ = right on the momentum space expression (fTHjl . This would lead to a certain relation 
between different parameters of the model, in our case it is 64^ = 2 — 6. If one takes all 
residues equal (i.e. with ^ = 0), it changes to the standard value 6 = 2. It is interesting to 
notice that for the parameters we have actually chosen A^ defined by 



as 
47r3 



A^ 



\2 ^ 

' 47r2 



2 I 2 

-m + rriQ 



(19) 



is equal to ^A^ = -0.2 GeV^ which is still "tachyonic'' 
-(0.07 - 0.09) GeV^ from 0. 



to be compared with — A^ 
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For the sake of completeness let us also cite the standard OPE answer for the discussed 
correlator (for the physical value Nc = 3) E2] '■ 



TT 3-27 ^ 81 



R^'{t} = 1^ ^ V^y T - ^:;-as{qq} r°logr /i^ + ... (20) 



where r is Euclidean time coinciding with xe in our case. We use the standard values 

{{gC;,)') = O.SGeV^ ; \{qq)\ = (250MeV)=' (21) 

It is worth noticing^ that the expression (j^UI) contains a contact term (j'^ju) which is not 
directly seen in the momentum space. This term's contribution is proportional to the sixth 
power of T in ()2()|1 . on the other hand, there is no ag in front of it [22] • Numerically according 

to [22! 

{u^^ddYu)/{qq)^ = -^(0.90 ± 0.15) 

We take this circumstance into account by means of numerical redefining fif, in (j^UI) accord- 
ingly, while still keeping the expression in the form ()20p . In fact, our main concern is to 
compare (fT^ with the lattice, while we need (j^U)) mostly for illustrative purposes. 



On Fig. 1.3 we compare the exact expression |T4|l with its own short distance expansion 
(I17j) and also with the conventional OPE result (pn|) where the three power terms have been 
kept in both cases. The Figure 1.3 shows striking difference from the Figs. 1.1, 1.2. Sizeable 
deviations of the standard OPE from p4|) start as early as at a; ~ 0.4 Fm. We have also 
plotted the first resonance contribution the answer should converge to at large distances. 
The sum over resonances smoothly interpolates between small and large distance regions. 



2.2 Axial channel 

We are now to consider the axial channel. Some care is needed here since one has to extract 
the pion contribution. There are two Lorentz structures in this channel 

VJ^Jiq)=tJd''xVJ^J{x)exp{tqx) = {q,qu - g^uq')M<l') + Qi^QuMq') (22) 

Contracting both sides with the tensor q^q^, one gets (see Appendix A): 

Qf^QuVl^'Jiq) = q^U2iq^) = 2(m„ + md){uu + dd) + (m„ + nidf P^^^q^) (23) 

We have put the bar over V^^\q^) in order to stress that we keep only the pion pole 
contribution in the pseudoscalar correlator in ()23|) (compare with the whole tower taken into 
account in (1311)). This is consistent with the use of Gell-Mann - Oakes - Renner relation 



^The authors are indebted to Prof. K.G.Chetyrkin for discussion of this point. 



in the standard form and such assumption has approximately the same level of accuracy as 
this relation itself. So we have 

P^'\-Q') = T^^^ (24) 

where as usual Q^ = — g^. The pion residue is fixed by PCAC (see Appendix A) as 

y^2 _ 2/^m^ 

^- - (m. + m,)2 ^''^ 

where we take /^ = 93 MeV, m^ = 140 MeV, mu + md = 11 MeV, so that A^ = 0.05 GeV^ 

In principle, one can address in the same framework the question about corrections 
to (|3H|) (and correspondingly to (jSHI), (EHl)) due to higher (nonchiral) or multipion states. 
This analysis is beyond the scope of the present paper. Notice the overall plus sign in (^^ 
fixed by the negative sign of the quark condensate, see Appendix A. One gets the following 
expression for 112 (g^): 

q'Me) = 4^ (26) 

while for ni(q'^) we have an expression analogous to (fTT|) : 

where we take tuq = Ma^ = 1.26 GeV while the residues coincide with those for the vector 
channel given by (J12p . We do not separate the lowest state (oi resonance) as we did for the 
vector channel since in this correlator the role of such state is played by pion. Moreover, 
as it has already been noticed the asymptotic value for residues is natural to expect to be 
rather close to the one for the lowest state. 



Contraction of Lorentz indices in ()22j) leads to 

V^''\q') = g.uVil'Jiq) = -Sq'Mq') + q'^q') (28) 

We have therefore 

Fourier thansformation to the coordinate space is straightforward, the answer is given by 
the following expression: 

R^''\^) = --^ fV^'M^z) + ^pt{z) (30) 

The following notations have been introduced: 7 = 771.,^ /mo, f = fTr/mQ. Notice the different 
sign of pion contribution and that of the tower. The correlator becomes negative at large 
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distances. On Figs. 2.1 and 2.2 expression (jHUj) is compared with the lattice data from j3] 
and from jl]. The short distance expansion of ()30|) is given by 



R^''\z) 



12 



2,3 



J 7 



^/ + plog(7^f + ... 



7 



+ 



+ 



1 + 



66-6-6^1 / + 



3-27 
while the standard OPE series can be written as 



36-2-62 
3^28 



z^logz^ + ... 



R^^U 



,^«.(r) {{gG'^.u?) 4^1l7r3 ,,2 6, 2 2^ 

1 + ^ ^v r +——as{qq) rMogr/i„ + ... 



TT 



3-27 



(31) 



(32) 



where we again denote xe as r and the discussion we had after ()2Tjl is also valid here. The 
comparison of (pUjl . (j!?T|) and (jHSI) is presented on Fig. 2.3. One can see that the standard 
OPE (as well as the short distance expansion (j!TT|) ) deviates significantly from the lattice 
data at distances larger than 0.4 Fm, while the full spectral density ()30|) reproduces the 
latter smoothly. 



2.3 Pseudoscalar channel 



This channel contains the light pion which makes it different from the others. The cor- 
responding Hamiltonian analysis should take the chiral nature of pion into account. The 
resulting spectrum contains the lightest excitation (pion) which is massless in the chiral 
limit plus a tower of massive states with the masses experiencing chiral shifts, which de- 
crease with n (see jS3])- Therefore we proceed as before, separating the lowest state (pion) 
and fixing the residues of all other poles by the duality requirement. The correlator reads as 



p(p)(g) =i [ d^x {0\T j^P\x)j^P\yy\0) exp{tq{x-y)) 



(33) 



with j^^\x) = u{x)i'j^d{x) . We take the polarization operator in the following form: 



V^P\-Q' 



XI 



+ E 



rip) 



Q' + ml ^,Q' + ml 



(34) 



We neglect the chiral shifts for higher states in what follows and take the spectrum of the 
tower in the form: 

m^ = m\ + rn^in — 1) ; n = l,2, ... (35) 

where mi = 1.3 GeV as for 7r(1300) state. The residues a^^^ should depend on n in such a 
way that the leading term is linear: 



aiP) = a^P^n + d^^^ + O (n-^) 



(36) 
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,(p) 



N, 



iTC^ 



m 



which provides 



This form is dictated by the quark-hadron duahty, postulated throughout the paper. How- 
ever, duahty requirement allows to fix only the leading term, a^P\ but not the subleading 
one a^^^ which plays a role of a free parameter in our analysis. For the coefficient a^'^ we 
obtain 

AT 

(37) 
(38) 

(39) 

(40) 



r^p\-Q' 



87r2 



Q'log 






The full expression for the correlator's ratio in coordinate space has the form: 

^2 h'^ia + l] 



TC 



where we use the following notation: 

7 = — , ^ = ^ , 
mi iTLi 



25 



h hz 

-z^Pbiz) + —pbiz) 



on ttp 
3m4 



nil 



rriiXE 



Notice the relative positive sign of the first and the third terms in the rhs of (jH^ . The 
parameters 7, L and b are fixed by their physical values while the parameter a is free in our 
approach, on physical grounds it is reasonable to expect it to be of the order of unity. On 
Figure 3.1 and 3.2 we plot the function R^'p\z) together with the lattice data for the choice 
a = 0. The Figure 3.4 shows the curves (jHn|l for the choices a = —1, 0, 2. Varying a slightly 
affects intermediate distance behavior of the correlator but does not touch short and long 
distance asymptotics. 

As above, we consider the short distance expansions. For (jH^ it reads 



IT 



R^^\z)='-L^^ 



-z^ + '^-zHog{^zf + ... 



+ 



h{a 



23 



z^ + 



+ 



+ 



1 + 



^i'" 



6-6^ / + 



36-2-62 
3^28 



z^Xogz^ + ... 



(41) 



where dots stay for higher powers of z"^. In contrast to (fTTjl . (p?T|l all non- negative powers of 
z^ are present in ()4H) . The conventional OPE answer in this channel has the form (see, e.g. 



R^P\i 



I 



r — 



3-2 



7 



If one includes the quadratic power correction in the OPE expansion 
relation is to be 



— a,(gg)^r''logr^/i2 + ... (42) 

, the corresponding 



TT 



6(a + l) 



1 



(43) 



Since a is a free parameter in our approach we can make no reasonable estimate of A^ from 
PH|l. The other way around, fixing the l.h.s. of (jlSjl by (fTU|) corresponds to the value a = 0.4. 
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The Figure 3.3 shows the exact answer (jHUj) together with the short distance asymptotics 
(JHj) and (jl21)- Analogously to the other channels, short- distance expansions break down at 
X ~ 0.4 Fm. 



2.4 Scalar channel 

This channel is presumably the most complex one since it requires some understanding of 
the ao meson nature,® in particular, the complicated dynamical problem of the role played 
by the four quark state admixture in these states. However, we do not address this set of 
questions here since the lattice data of we are going to compare with are obtained in 
quenched approximation. We also take the set of data corresponding to the largest value of 
current quark mass used in simulations ^. We also have in mind large Nc limit, where four 
quark admixture vanishes. We write therefore 

'P'''i-Q') = to^ (44) 

n=0 V "T "'n 

with the corresponding quark-hadron duality condition 

p(^)(-g^)=^^°^^gMog(^) (45) 

The expression for residues coincides with fl3(ij) : 

flW = a'^'^n + a^'^ + C ( r ) (46) 

and 

a(^) = 



vith 




-,Q'ic 


'g 




a^'^n + a(^) 


+ 




) 


a(^) 


N, 


m^ 







(47) 



The result for the ratio R^'^\z) is 



R^^\z) = ^-^^^z^p,iz) + ^-^Mz) (48) 



The notation is standard: 



Stt a, , m 



2 



3m'^ nin 



rriQXE 



Notice the difference between (HHjl and (jHOI)- In the latter case the vr - meson pole stay away 
from the tower, but the numbering of states still respects the fact that the pion is the lowest, 
n = state. This is in line with j3]. In the scalar channel the tower itself starts from the 



®Let us remind that we work with 7=1 case only, so we have no mixing with the pure glue states and 
do not consider /q meson. 
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lowest state. This causes 6^ (a + 1) — 6 factor in (jH^ instead of b{ab — 1) in PHj). Also the 
parameter b is defined differently. One may check that if mass and residue of the lowest state 
are chosen in such a way that this state also belongs to the equidistant tower, the expression 
(I48|) identically coincides with (|39|) . as it should be. 

The standard OPE answer is given by j2] 

R^'\r) = 1 + M^!! r^ + ^a.(gg)^ r« logr^^ + •.• (49) 



while the short distance expansion of (pSj) takes the form 

ba — 1 



R^'h 



23 



z' + 



+ 



+ 



1 + ^(66-6-6^) z'+ ^\^^,^\ 'logz' + 



(50) 



The Figure 4.1 shows the curve defined by (J4H|) for different values of a against the 
data from [3]. Unfortunately, the data in this channel are not reliable enough, which can 
be clearly seen on the Figure 4.2, where the lattice correlator from j^ becomes negative at 
X ~ 0.5 Fm due to lattice artefacts, in contradiction with positivity of the spectral density.''' 
The short distance expansions explode even earlier (see Figure 4.3). 



3 Discussion of the results and conclusion 

We find the results of comparison of the lattice data and the model to be quite remarkable. 
It is worth stressing that we have worked with the model of two relativistic quarks connected 
by the string with the tension a, which is the only dimensionful parameter of the model. 
Using this picture, first, the lowest resonance masses were computed in all channels J2] 
and found to be in reasonable agreement with their experimental values. Second, the quasi- 
classical asymptotic for the mass spectrum was calculated J7] which has the same pattern 
for all channels. The latter must provide exact quark-hadron duality, which was explicitly 
shown in the vector case in ^^1 ISO] • We have plugged these two ingredients into the cor- 
responding spectral density and compared the results with the available lattice simulations 
13 H] . Another comparison is made for the standard S VZ short distance OPE expansion 
and the corresponding expansion provided by our spectral density. It is shown that they 
strongly deviate from the full curve (and from the lattice data) for the distances larger than 
0.35 — 0.45 Fm, as one could expect. 

The main lesson is that one needs quite a few inputs (correct lowest resonance mass, 
correct lowest resonance residue and correct asymptotic behavior dictated by quark-hadron 

^The authors are grateful to Prof. De Grand for discussion of the problems in this channel. 
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duality) to reproduce lattice data in a reasonable way. We have taken these inputs from 
the large Nc model of QCD string with quarks at the ends [TTj. The absence of precise 
quantitative agreement between our results and the lattice should not disappoint since many 
effects have been ignored (1/iVc effects, perturbative exchanges etc), notice also that the 
two sets of lattice data we have used do not agree with each other on quantitative level. 
On the other hand, qualitative agreement is rather good and certainly better than that of 
the standard short distance OPE in coordinate space. Notice that we have not performed 
any fitting of the lattice data, the latter was compared with the curves (|T4|l . pO|l . p9|) . (J48|l 
computed independently. If we had fitted the data with these expressions, the agreement 
would have been much better. 

The next logical step is to compare the predictions of large Nc QCD string model with 
the data taken from real experiments, i.e. the data on r - decay rates obtained by ALEPH 
and OPAL collaborations. The accuracy of the data is much higher than what is possible to 
reach on the lattice, and this fact provides a challenging task for theory. The r - decay results 
have been already addressed in the framework of conventional OPE [SI] and in the instanton 
model J2S1- The issues of large Nc limit interpretation of the data have also been discussed 

We are going to analyse these data in our framework in a separate publication. 
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Appendix A 

This Appendix contains a textbook material (see, e.g. ^), however we find it useful for the 
sake of completeness to rederive here the equation (|23|) in our conventions. The following 
expressions for charged axial and pseudoscalar currents have been introduced: 

jI^\x) = u{x)'jf,'j5d{x) ; j^P\x) = u{x)i'j5d{x) (51) 

Using Dirac equation 

(z7^(9^ - igA^ix)) - m)^ix) = (52) 

one gets the well known expression for PCAC: 

dJ;\x) = {m^ + m,)j^P\x) (53) 
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Next, one defines the following matrix elements of the currents between vacuum and one 
pion state 

{0\jl'\x)\n{p)) = tV2 Up^exp{-ipx) ; {0\j^^\x)\7r{p)) = X^exp{-ipx) (54) 

Our normalization coincides with P] but \/2 different from [S]. Correspondingly, we have 
/^ = 93 MeV instead of 131 MeV value for /^ used in jS]. 



Combining (jH!?jl and (jH^ one gets 



V^fnin: 



A. = ^"^7"'' (55) 

Consider now the current-current correlator in momentum space 

Vl^J{q)=zJd'x{0\Tjl^\x)£^\yy\0) exp{tq{x - y)) (56) 

Making use of free equal-time anticommutation relations for the quark fields, it is straight- 
forward to obtain: 

^quV^^iq) = {"^u + mrf)2p(p)(g2) + 2(m„ + mi)l,uu + dd) (57) 

where V^^\q) is given by (jHHj) . Assuming that in the chiral limit one pion state contribution 
dominates in the pseudoscalar correlator at g^ -^ 0: 

ml 

one gets 

2ml i^u + dd) + (m„ + md)\l = (59) 

which together with (|55|) is nothing but the Gell-Mann - Oakes - Renner relation 

(m„ + ma) {uu + dd) = -f^ml (60) 

up to the terms of higher order in m^. Taking into account the definition (|22p we come to 
the expression (PH|) used in the main text. 



Appendix B 

We collect here the basic formulas used in the main text. The standard metric conventions, 
correlator definitions, Fourier transforms are adopted. The free propagator is given by 

- f d^q q + m ( d \ m? 

m^)m) = V (^ ^^^:^2 exp(-.gx) = (^^7.^ + "^j 4^^ ^^ (^^) (61) 
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for space-like x, where z = m\/ —x^. Four vectors in Minkowskii and Euclidean spaces are 
defined as follows: 

1 = (?o, q) ; q'^ = qli = ql-^ ; Q^ = ql = -qli (62) 

and we often omit the subscripts ii^, M if it is not misleading. Notice that the Wick rotation 
changes the free propagator sign : 

The Fourier transform for the free propagator is given by (see, e.g. i2I])- 

/•I TT^ 

J ^^^ T^ exp(igx) ^ ^y (-?^) log(-?^) (64) 



Let us calculate the function P{x) we have made use of in the main text. It is defined 
as follows: 

Pi.^) = 3^ / jTT^ Q^ Y. — —, 2 exp(-igx) (65) 

J (27r)4 ^^-q2 + rnl 

where m^ = rriQ + rn^n. The sum over n is divergent and the divergent (g^ - independent) 
constant term is usually eliminated by renormalization n(g^) -^ n(g^) — 11(0) in momentum 
space. We proceed in a different way here, simply interchanging the summation over n and 
integration over q for x ^ 0, having in mind condition (jU)) then a; ^ 0. We thus systematically 
disregard all terms proportional to 6{x) or derivatives of the delta function. In this way one 
obtains 

^W-3(-^)Ej^A-.(,„,.^) (66) 



We switch to the Euclidean notations in the summand (jUUI) . with xe = \f—x^ > 0. Using 
integral representation for the McDonald function 



oo ,. / 2^ 



^ ^ 



and performing the summation over n 



x|m^ \ ^^P V 4t 



4t 
we obtain 



^exp -i£^ = V^-f- (68) 



n=o \ /I — exp I — -. 



2 9 

4t 



^w-dbO/f-p(-"4^) '•''' 
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where 



E 



x%ml' 



^°l"j^?o'""'''^ « 



exp 



4t 



1 — exp 



x^^m? 



it 



ml + (m^ — r?7,QJ exp 



(70) 



To represent the results in more convenient form we use the notation z = xeitiq for 
dimensionless distance and b = rn^ /rn^ for the quanta/gap ratio. Changing integration 
variables t -^ u such that -iut = z^, one gets 



P{x) 



du 



92 



^ 



4 

^r^xp 
x^ 



x^rr^ 
Au 



exp(— Mj 



'1 — exp(— 6m))^ 



(l + exp(-6M)(6-l)) 



Performing differentiation, we finally obtain 

Pix) = - 



3mQ 



167r2 



Pb[z) 



(71) 
(72) 



where the universal function pb{z) has the following form 



du f z^ \ l + exp(-6M)(6- 11 
M^J=/^exp(^---.j 

^ ^ 



(1 — exp(— 6n))^ 

This result has been used in the main text. 

For scalar and pseudoscalar channels one also needs the function P[x) defined as 

d^q ^ 1 



(73) 



P{x) = -i 



W Sv + K 



exp(— iga;) 



(74) 



where again m'^ = rn^ + m?n. Performing completely analogous calculations one gets 



l + exp(-6M)(6- 1) 
(1 — exp(— 6m))2 



pb{z) = Ait x^ P{x) = j du exp — u 



(75) 



The short-distance asymptotics oi pb{z) can be easily found from the definition (fTSj). it 
is given by: 



/ X 128 1 / 2 b\ 1 36-2-62 



(76) 



where the dots denote constant (z^ -independent) term and positive powers of z"^. 
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Figure captions 

Fig. 1.1 Lattice data from [T for the vector channel correlator vs expression (14). 

Fig. 1.2 Lattice data from [4. for the vector channel correlator vs expression (14). 

Fig. 1.3 Expression (14) (solid curve) together with its own short distance expansion (17) 
(dashed-dotted curve) and short distance OPE expansion (20) (dashed curve). The 
lightest resonance contribution is also shown (double dotted curve). 

Fig. 2.1 Lattice data from [^ for the axial channel correlator vs expression (30). 

Fig. 2.2 Lattice data from [4^ for the axial channel correlator vs expression (30). 

Fig. 2.3 Expression (30) (solid curve) together with its own short distance expansion (31) 
(dashed-dotted curve) and short distance OPE expansion (32) (dashed curve). 

Fig. 3.1 Lattice data from [3 for the pseudoscalar channel correlator vs expression (39). 

Fig. 3.2 Lattice data from [4 for the pseudoscalar channel correlator vs expression (39). 

Fig. 3.3 Expression (39) (solid curve) together with its own short distance expansion (41) 
(dashed-dotted curve) and short distance OPE expansion (42) (dashed curve). 

Fig. 3.4 Expression (39) for different values of the parameter a: a = 2 (upper curve), a = 
(solid curve) and a = — 1 (lower curve). 

Fig. 4.1 Lattice data from [HI for the scalar channel correlator vs expression (48) for different 
values of a: a = 2 (upper curve), a = 1 (middle curve) and a = (lower curve). 

Fig. 4.2 Lattice data from [T for the scalar channel correlator vs expression (48). 

Fig. 4.3 Expression (48) (solid curve) together with its own short distance expansion (50) 
(dashed-dotted curve) and short distance OPE expansion (49) (dashed curve). 
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